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Abstract 

Results from a recent analysis of the zero-skewness generalized parton distribu¬ 
tions (GPDs) for valence quarks are reviewed. The analysis bases on a physically 
motivated parameterization of the GPDs with a few free parameters adjusted to the 
nucleon form factor data. The Fourier transforms of the GPDs representing quark 
densities in the impact parameter plane, as well as moments of the GPDs are also 
discnssed. The 1/x moments in particular form the soft physics input to Compton 
scattering off protons within the handbag approach. The Compton cross section 
evaluated from this information is found to be in good agreement with experiment. 

Talk presented at the XI workshop on High Energy Spin Physics (SPIN 2005), 
JINR Dubna, September 2005 


1 Introduction 

In recent years we have learned how to deal with hard exclusive reactions within QCD. In 
analogy to hard inclusive reactions the process amplitudes factorize in partonic subpro¬ 
cess amplitudes, calculable in perturbation theory, and in GPDs which parameterize soft 
hadronic matrix elements. In some cases rigorous proofs of factorization exist. For other 
processes factorization is shown to hold in certain limits, under certain assumptions or is 
just a hypothesis. The GPDs which are dehned by Fourier transforms of bilocal proton 
matrix elements of quark held operators [1], describe the emission and reabsorption of 
partons by the proton. For equal helicities of the emitted and reabsorbed parton the 
structure of the nucleon is described by four GPDs, termed H, H, E and E, for each 
quark havor and the gluons. The GPDs are functions of three variables, the momentum 
transfer from the initial to the hnal proton, t, a momentum fraction x and the skewness. 
The latter variable is related to the difference of the individual momentum fractions 
the emitted and reabsorbed partons carry. The GPDs are subject to evolution and, hence, 
depend on the factorization scale p, too. They satisfy the reduction formulas 

^ = 0, t = 0, p) = q{x, p), H'^{x, ^ = 0, f = 0, /i) = Aq{x, p), (1) 

i.e. in the forward limit H and H reduce to the usual unpolarized and polarized parton 
distributions (PDFs), respectively. Another property of the GPDs is the polynomiality 
which comes about as a consequence of Lorentz covariance 
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where [m/2] denotes the largest integer smaller than or equal to m/2. Eq. (2) holds 
analogously for the other GPDs and, for m = 1, implies sum rules for the form factors of 
the nucleon, e.g. 

Fi{t) = h%{t) = dxH\x,^,t,fx). (3) 

Reinterpreting as usual a parton with a negative momentum fraction x as an antiparton 
with positive x {H^{x) = —H'^{—x)), one becomes aware that only the difference of the 
contributions from quarks and antiquarks of given flavor contribute to the sum rules. 
Introducing the combination 






( 4 ) 


for positive x which, in the forward limit, reduces to the usual valence quark density 
qv{x) = q{x) — q{x), one hnds for the Dirac form factor the representation 

= eu{d) dxH^{x,^,t,iJ,) + ed(u)J^ dx H^{x,^,t, n). (5) 

Here, is the charge of the quark q in units of the positron charge. Possible contributions 
from other flavors, s — s or c — c, are neglected in the sum rule (5). A representation 
analogous to (5) holds for the Pauli form factor with H being replaced by E. 

The isovector axial-vector form factor satishes the sum rule {H^{x) = H^{—x)) 


i(i) = C, t: h) - h) + 2 [h^{x, t, /r) - H‘^{x, t, /i)] } , (6) 


At least for small t the magnitude of the second integral in Eq. ( 6 ) reflects the size of 
the flavor non-singlet combination Au{x) — Ad{x) of forward densities. This difference 
is poorly known, and at present there is no experimental evidence that it might be large 
[2]. For instance, in the analysis of the polarized PDFs performed by Bliimlein and 
Bottcher Ref. [3] it is zero. In a perhaps crude approximation the second term in ( 6 ) 
can be neglected. These simplihcations to the sum rules (5) and ( 6 ) do not imply that 
the nucleon is assumed to consist solely of valence quarks. Sea quarks are there but the 
virtual photon that probes the quark content of the nucleon, sees only the differences 
between quark and antiquark distributions which are likely small. 


2 A determination of the zero-skewness GPDs 

Since the GPDs cannot be calculated from QGD with a sufficient degree of accuracy at 
present we have either to rely on models or to extract them from experiment as it has 
been done for the PDFs, see for instance Refs. [3, 4]. A Erst attempt to extract the 
GPDs phenomenologically has been carried through in Ref. [5] (for a similar analysis 
see Ref. [ 6 ]). The idea is to exploit the sum rules (5) and ( 6 ) and to determine the 
GPDs from the data [7] on the nucleon form factors, Fi and F 2 for proton and neutron 
as well as i/ 4 . Since the sum rules represent only the hrst moments of the GPDs this 
task is an ill-posed problem in a strict mathematical sense. Inhnitely many moments 
are needed to deduce the integrand, i.e. the GPDs, unambiguously from an integral. 
However, from phenomenological experience with particle physics one expects the GPDs 
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to be smooth functions and a small number of moments will likely suffice to fix the 
GPDs. The extreme - and at present the only feasible - point of view is that the lowest 
moment of a GPD alone is sufficient [5]. Indeed, using recent results on PDFs [3, 4] 
and form factor data [7] in combination with physically motivated parameterizations of 
the GPDs, one can carry through this analysis. Needless to say that this method while 
phenomenologically successful as will be discussed below, does not lead to unique results. 
Other parameterizations which may imply different physics, cannot be excluded at the 
present stage of the art. 

The sum rules (5), (6) are valid at all ^ but guessing a plausible parameterization of 
a function of three variables is a nearly hopeless task. Ghoosing the special value = 0 
for which the emitted and reabsorbed partons carry the same momentum fractions, has 
many advantages. One exclusively works in the so-called DGLAP region where ^ < x. In 
this region parton ideas apply and the Fourier transform of a GPD with respect to the 
momentum transfer A (A^ = —t) has a density interpretation in the impact parameter 
plane [8]. Moreover, as shown in Refs. [9, 10], wide-angle exclusive processes are controlled 
by generalized form factors that represent 1/x moments of zero-skewness GPDs. The 
parameterization that is exploited in Ref. [5] combines the usual PDFs with an exponential 
t dependence (the arguments = 0 and fx are omitted in the following) 

= qy{x)exp[tfg{x)], (7) 

where the profile function reads 

fq{x) = [a'\og{l/x) + Bq]{l - + Aqx{l - x)"^. (8) 

This ansatz is motivated by the expected Regge behavior at low —t and low x [11] {a' is 
the Regge slope for which the value 0.9 GeV^ is imposed). For large —t and large x, on 
the other hand, one expects a behavior like fq^l — x from overlaps of light-cone wave- 
functions [9, 12]. The ansatz (7), (8) interpolates between the two limits smoothly ^ and 
allows for a stronger suppression of fq in the limit x ^ 1. It matches the following criteria 
for a good parameterization: simplicity, consistency with theoretical and phenomenologi¬ 
cal constraints, stability with respect to variations of PDFs and stability under evolution 
(scale dependence of the GPDs can be absorbed into parameters). 

Using the GTEQ PDFs [4], one obtains a reasonable fit of the ansatz (7), (8) to the 
data on the Dirac form factor (from —f = 0 np to ~ 30 GeV^) with the parameters 

Bu = Bd = (0.59 ± 0.03) GeV^ , 

Au = (1.22 ± 0.020) GeV^ , Ad = (2.59 ± 0.29) GeV^ , (9) 

quoted for the case n = 2 and at a scale of /i = 2 GeV. In Fig. 1 the resnlts for 

H!^ are shown at three valnes of t. While at small —t the behavior of the GPD still 
reflects that of the parton densities it exhibits a prononnced maximnm at larger valnes of 
—t. The maximnm becomes more prononnced with increasing —t and its position moves 
towards x = 1. In other words only a limited range of x contribntes to the form factor 

snbstantially. The qnality of the fit is very similar in both the cases, n = 1 and 2; the 

resnlts for the GPDs agree well with each other. Snbstantial differences between the two 

^The parameter Bq is not needed if a' is freed. A fit to the data of about the same quality and with 
practically the same results for the GPDs is obtained for a' ~ 1.4. 
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Figure 1: Results for H^{x,t) at the scale /i = 2 GeV and for u = 2 obtained in [5]. In 
each of the regions indicated by dashed lines, 5% of the total value of the sum rule (5) 
are accumulated. 

results only occur for very low and very large values of x, i.e. in the regions which are 
nearly insensitive to the present form factor data. It is the physical interpretation of the 
results which favours the fit with n = 2. Indeed the average distance between the struck 
quark and the cluster of spectators becomes unphysical large for x —>■ 1 in the case n = 1; 
it grows like ~ (1 — x)“^ while, for n = 2, it tends towards a constant value of about 0.5 
fm [5]. 

The analysis of the axial and Pauli form factors, with parameterizations analogous to 
Eqs. (7), (8), provides the GPDs H and E. In general they behave similar to H. Note¬ 
worthy differences are the opposite signs of (£’“) and H!^ {Ef) and the approximately 
same magnitude of E^ and E^ at least for not too large values of —t. For H!^ and 
on the other hand, the d-quark contributions are substantially smaller in magnitude than 
the M-quark ones, see Fig. 1. Since there is no data available for the pseudoscalar form 
factor of the nucleon the GPD E cannot be determined this way. 


3 Moments of the GPDs 

Having the GPDs at disposal one can evaluate various moments, some of them are dis¬ 
played in Fig. 2. Gomparison with recent results from lattice QGD [13] reveals remarkable 
agreement of their t dependencies given the uncertainties of the GPD analysis [5] and of 
the lattice calculations. It is to be stressed that the results presented in Refs. [13] are pro¬ 
vided for a scenario in which the pion has a mass of 600 to 800 MeV. The extrapolation 
to the chiral limit is therefore problematic and may entail failures of the normalization 
of the moments. An interesting property of the moments is that the u and d quark con¬ 
tributions decrease with different rates at large —t. Because at large —t the dominant 
contribution to the form factors comes from a narrow region of large x, see Fig. 1, one 
can use the large x approximations 

qyr^{l-xf\ /g~Ag(l-x)^ (10) 
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Figure 2: Left: The first three moments of scaled by t^. The error bands denote the 
parametric uncertainties. Right: The lowest moment of scaled by The shaded 
band represents the results obtained in Ref. [5]. Data are taken from [7, 15]. 


and evaluate the sum rule (5) in the saddle point approximation. This leads to 

( \ -i/n 

, ( 11 ) 

where Xg is the position of the saddle point. For n = 2 the saddle point lies within the 
region where the GPD is large ^ and using the the CTEQ values for jSq ~ 3.4 and 
f3d — 5 [4]) one obtains a drop of the form factor slightly faster than t~‘^ while the d- 
quark form factor falls as t~^. Strengthened by the charge factor the u-quark contribution 
dominates the proton’s Dirac form factor ^ for —t larger than about 5 GeV^, the d-quark 
contribution amounts to less than 10%. This is what can be seen in Figs. 2. The power 
behavior bears resemblance to the Drell-Yan relation [14]. In fact the common underlying 
dynamics is the Feynman mechanism The Drell-Yan relation is, however, an asymptotic 
result (x —> 1, f ^ —cxo) which bases on the assumption of valence Fock state dominance, 
i.e. on the absence of sea quarks. The different powers for u and d quarks signal that the 
asymptotic region where the dimensional counting rules apply, has not yet been reached. 

One may object to the different powers that they are merely a consequence of the 
chosen parameterization (7), (8). However, this is likely not the case. The moments 
can directly be extracted from the form factor data [7]. The experimental results for the 
d-quark moment are shown in Fig. 2. The sharp drop of it for increasing —t is clearly 
visible. The small deviations between data and the moment obtained in [5] are caused 
by the very precise but still preliminary GLAS data [15] for These data have been 

^For n = 1 this is only the case for —t larger than 30 GeV^. The onset of the corresponding power 
behaviour of the form factors occurs for larger —t. 

•^This implies the ratio F^/Ff = e^/e^ at large —t. 

^The Feynman mechanism applies in the soft region where 1 — x ~ A/-%^ and the virtualities of the 
active partons are ~ A-%^ (A is a typical hadronic scale). 
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Figure 3: Left: The electric form factor of the neutron versus t [16]. The lines represent 
different extrapolations, see text. Right: Sketch of flavor segregation in a transversally 
polarized proton. 


utilized in the extraction of the moments hf q but not in the analysis in Ref. [5]. Worst 
measured of the four form factors is the electric one of the neutron. The recent data [16] 
on it are shown in Fig. 3. Above 1.5 GeV^ no data is available as yet and extrapolations to 
larger —t are to be used. The solid line has been used in Ref. [5] and in the determination 
of hf 0 from data. If the electric form factor is smaller than that extrapolation (e.g. the 
dashed line) falls even faster. Thus, except is much larger than expected, the 
rapid decrease of hfg seems to be an experimental fact. Future JLab data on will 
settle this question. 

A combination of the second moments of H and i? at t = 0 is Ji’s sum rule [17] 
which allows for an evaluation of the valence quark contribution to the orbital angular 
momentum the quarks inside the proton carry 


m) 


1 

2 


dx 


xEl{x, t = 0) + xq^{x) - Aqy{x) 


The analysis of the GPDs leads to 


{L“) = -(0.24-0.27), {Lfi = 0.15-0.19, 


( 12 ) 


(13) 


for the valence quark contributions to the orbital angular momentum at a scale of /i = 
2 GeV. 

While the parton distributions only provide information on the longitudinal distri¬ 
bution of quarks inside the nucleon, GPDs also give access to the transverse position 
distributions of partons within the proton. Thus, the Fourier transform of H 

r 

= J = (14) 

gives the probability of hnding a valence quark with longitudinal momentum fraction x 
and impact parameter b = [b^, b^) as seen in a frame in which the proton moves rapidly 
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in the Z direction. Together with the analogous Fourier transform of E'^{x,t) one can 
form the combination [rUp being the mass of the proton) 

d 

q^{x, b) = qy{x, b)-— el{x, b) , (15) 

rrip ob^ 

which gives the probability to hnd an unpolarized valence quark with momentum fraction 
X and impact parameter b in a proton that moves rapidly along the Z direction and is 
polarized along the X direction [8]. As shown in Ref. [5], for small x one observes a very 
broad distribution while at large x it becomes more focused on the center of momentum 
dehned by J2iXihi = 0 (J2iXi = !)• In a proton that is polarized in the X direction 
the symmetry around the Z axis is lost and the center of the density is shifted in the 
Y direction away from the center of momentum, downward for d quarks and upward for 
u ones. In other words, a polarization of the proton induces a flavor segregation in the 
direction orthogonal to the direction of the polarization and the proton momentum, see 
Fig. 3. This effect may be responsible for certain asymmetries as, for instance, that one 
observed in p | p — tt^X. 


4 Wide-angle scattering 


As mentioned previously the analysis of the GPDs gives insight in the transverse distri¬ 
bution of quarks inside the proton. However, there is more in it. With the = 0 GPDs 
at hand one can predict hard wide-angle exclusive reactions like Gompton scattering off 
protons or meson photo- and electroproduction. For these reactions one can work in a 
so-called symmetrical frame where the skewness is zero. It has been argued [9, 10] that, 

for large Mandelstam variables (s, —t, —u 3> A^), the amplitudes for these reactions 

factorize into a hard partonic subprocess, e.g. Gompton scattering off quarks, and in form 
factors representing l/x-moments of zero-skewness GPDS (see Fig. 4). For Gompton 
scattering these form factors read 

Rvit) ~ Y. 4 f > 

, ^ JO X 
q=u,d 

RaJ) ~ Y 4 [ 

. . ^ -'0 X 

q=u,d 

_ rl fix 

RtJ) - Y4 —Rv{x,t)- (16) 

q=u^d 


In these expressions contributions from sea quarks have been neglected An analogous 
pseudoscalar form factor related to the GPD E, decouples in the symmetric frame. Nu¬ 
merical results for the Gompton form factors are shown in Fig. 4. To leading-order of 
pQGD the Gompton cross section reads[9, 18] 


da 

dt 


dd 

dt 


1 r 

2 - 


R^Jt) + '^:^RtJ) + R\{t) 


®An estimate of the sea quark contributions may be obtained by using the ansatz (7) with the same 
profile function for the sea quarks as for the valence ones but replacing the valence quark density with 
the CTEQ [4] antiquark ones. The so estimated contributions are indeed small. 
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Figure 4: Left: The Compton form factors, scaled by evaluated from the GPDs deter¬ 
mined in Ref. [5]. The bands represent the parametric uncertainties of the form factors. 
Right: The Compton cross section versus the c.m.s. scattering angle 6 ^ at s = 6.9, 8.9, 11 
and 20 GeV^. The error bands are explained in the text. Data are taken from Ref. [19]. 


us 







(17) 


where da/dt is the Klein-Nishina cross section for Compton scattering off massless, point- 
like spin-1/2 particles of charge unity. Inserting the Compton form factors (16) into Eqs. 
(17), one can predict the Compton cross section in the wide-angle region. The results 
for sample values of s are shown in Fig. 4 and compared to recent measurements from 
JLab [19]. The inner bands of the predictions for da/dt reflect the parametric errors of 
the form factors. The outer bands indicate the size of target mass corrections, see Ref. 
[20]. In order to comply with the kinematical requirements for handbag factorization, at 
least in a minimal fashion, predictions are only shown for —t and —u larger than about 
2.5 GeV^. Fair agreement between theory and experiment is to be seen. Next-to-leading 
order QGD corrections to the subprocess have been calculated in Ref. [18]. They are not 
displayed in (17) but taken into account in the numerical results. 

The handbag approach also applies to wide-angle photo- and electroproduction of 
pseudoscalar and vector mesons [21]. The amplitudes again factorize into a parton-level 
subprocess, yg —>■ Mq now, and form factors which represent l/x-moments of GPDs. 
Their flavor composition reflects the valence quark content of the produced meson. There¬ 
fore, the form factors which occur in photo-and electroproduction of pions and p mesons, 
can also be evaluated from the GPDs given in Ref. [5]. 

One may also consider the time-like process 77 •«-*• pp within the handbag approach 
[22]. Similar form factors as in the space-like region occur but they are now functions of 
s and represent moments of the pp distribution amplitudes, time-like versions of GPDs. 
With sufficient data on the time-like electromagnetic form factors at disposal one may 
attempt a determination of the time-like GPDs. 













5 Outlook 


Results from a first analysis of GPDs at zero skewness have been reviewed. This analysis, 
has been performed in analogy to those of the usual parton distributions, rests upon a 
physically motivated parameterization of the GPDs htted to the availablejiucleon form 
factor data. The analysis provides results on the valence-quark GPDs H, H and E. The 
distribution of the quarks in the impact parameter plane transverse to the direction of 
the nucleon’s momentum can be evaluated from them. Polarizing the nucleon induces a 
flavor segregation in the direction orthogonal to the those of the nucleon’s momentum and 
of its polarization. The average orbital angular momentum of the valence quarks can be 
estimated from the obtained GPDs, too. Within the handbag approach the soft physics 
input to hard wide-angle exclusive reactions is encoded in specihc form factors which 
represent 1/x moments of zero-skewness GPDs. Using the GPDs determined in Ref. [5], 
one can evaluate these form factors and, for instance, predict the wide-angle Gompton 
cross section. The results are found to be in good agreement with experiment. 

An analysis as that one performed in Ref. [5] can only be considered as a first attempt 
towards a determination of the GPDs. It needs improvements in various aspects. High 
quality data on the form factors at larger t are required in order to stabilize the parame¬ 
terization. JLab will provide such data in the near future. As already mentioned GLAS 
[15] will come up with data on up to about 5 GeV^. will be measured up to 
about 3.3 GeV^ next year and up to 9 GeV^ in 2007. The upgraded Jlab will allow 
measurements of the nucleon form factors up to about 13 GeV^. Data on the axial form 
factor are needed in order to improve our knowledge of H. 

Up to now only the sum rules (5), (6) have been utilized in the GPD analysis. As 
mentioned above, in this situation parameterizations of the GPDs are required with the 
consequence of non-unique results. An alternative ansatz [5] is for instance 




qv[x) 


1 - tfq{x)/p 



(18) 


Although reasonable fits to the form factors are obtained with (18) for p > 2.5, it is physi¬ 
cally less appealing than the parameterization (7): the combination of Regge behavior at 
small X and t with the dynamics of the Feynman mechanism is lost. The resulting GPDs 
have a broader shape and H{x = 0,t) remains hnite. Thus, small x also contribute to 
the large-t form factors. Higher moments are needed in order to lessen the dependence 
on the chosen parameterization for the GPDs. Such moments can be provided by lattice 
gauge theories. The present lattice results [13] are however calculated in scenarios where 
the pion is heavy (typically 600 — 800 MeV) and the extrapolation to the chiral limit is 
uncertain. Obviously such results are inappropriate for use in a GPD analysis. In a few 
years the quality of the moments from lattice gauge theories may suffice. It is also tempt¬ 
ing to use the data on wide-angle Gompton scattering [19] which imply information on 
1/x moments of the GPDs. However, even at the highest measured energy, s = 11 GeV^, 
the Gompton cross section may still be contaminated by power corrections rendering its 
use in a GPD analysis dubious. For data at, say, s = 20 GeV^ the situation may be 
different. Instead of a parameterization one may use the maximum entropy method [23] 
for the extraction of the GPDs from the sum rules. This method has not yet been utilized 
for this purpose. 

Up to now only the zero-skewness GPDs have been determined. Although the sum 
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Figure 5: Experimental information on GPDs. 


rules (5), (6) are valid at all values of the skewness it seems a hopeless task to extract 
functions of three variables from them. Additional information is demanded in the case of 
7 ^ 0 and will be provided by deeply virtual exclusive scattering (DVES). The skewness is 
related to Bjorken-x in DVES and, hence, non-zero. At large but small t the amplitudes 
for DVES are given by a convolution of GPDs and appropriate subprocess amplitudes. 
The convolution possesses poles at x = ±^. For Gompton scattering the pole terms can 
directly be isolated in the interference region with the Bethe-Heitler process and will lead 
to a nearly model-independent determination of For x ^ ^ parameterizations 

of the GPDs will most likely be needed. At present a number of experiments on DVES are 
running (HERA, GOMPASS, HERMES, JLab). In Fig. 5 kinematical regions are shown 
as shaded areas in which information on GPDs is currently accumulated. A lot of work 
is ahead of us before we can say that we have a fair knowledge on the GPDs. 
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